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THE ISOTOPE EFFECT IN (5-ZIRCONIUM 
by Donald Graham and Earl R. Hanes, Jr. 


Lewis Research Center 
SUMMARY 

89 95 

The self-diffusion isotope effect has been determined for Zr and Zr diffusing in 
|3-zirconium. The migration rate is independent of the tracer mass within experimental 
error at 942°, 1186°, 1784°, and 1815° C. The maximum value of diffusion coefficient 
ratio Dgg/Dgg is 1. 002±0. 005. 

An anneal at 1750° C for 8 hours immediately followed by the deposition and diffu- 
sion of niobium-95 in the /3-phase at 1025° C resulted in a threefold reduction of the dif- 
fusion coefficient over that measured conventionally by depositing the tracer at room 
temperature in the a -phase. 

It is tentatively suggested that the anomalous diffusion characteristics of zirconium, 
titanium, and hafnium can be accounted for by the presence of a dense network of short 
circuit paths introduced into the j3- phase during the a/fi transformation. These de- 
fects appear to be relatively stable and can only be reduced in concentration by a pro- 
longed anneal close to the melting temperature. 


INTRODUCTION 

Self- diffusion in the "normal” body-centered-cubic metals such as tantalum (ref. 1) 
and molybdenum (ref. 2) may be represented by an Arrhenius or thermal activation re- 
lation D t = Dge _< ^/RT which describes the temperature dependence of the diffusion 
coefficient in terms of a temperature independent energy of activation Q which is re- 
lated empirically to the melting temperature (ref. 3). The frequency factor Dq takes 
values from 0. 1 to 10 square centimeters per second (ref. 4) and can be adequately de- 
scribed in terms of the thermodynamic and geometric properties of the lattice. 

In the anomalous metals titanium (ref. 5), zirconium (ref. 6), and hafnium (ref. 7), 
an Arrhenius representation gives a temperature dependent value of Q in the /3-phase 
which increases from approximately 1. 30 eV just above the transition from close packed 


hexagonal to body centered cubic to 3. 00 eV at the melting point compared to a Q of 
3. 00 to 3. 50 eV suggested by the melting rule. In the temperature range 800° to 1400° C 
the frequency factor is at least an order of magnitude smaller than that for a normal 
metal and is characteristic of dislocation or grain boundary migration rather than 
volume diffusion (ref. 8). 

Graham suggested in an earlier paper (ref. 9) that two competitive mechanisms 
could explain the anomalous temperature dependence of diffusion. The approach towards 
normal values of Q and Dq at high temperature indicated that perhaps a single vacancy 
exchange mechanism could only effectively dominate the diffusion behaviour at high tem- 
perature because of the presence of a low Q mechanism which existed throughout the 
/3-phase. 

The isotope effect measurement provides a test for this hypothesis. The recent 
measurement of the isotope effect for self-diffusion in body centered cubic iron (ref. 10) 
is consistent with diffusion via single vacancy exchange and is equal to 0. 72 within the 
experimental error. If single vacancy exchange occurs in /3-Zr at high temperature, we 
should observe an isotope effect equal to 0.72. At lower temperatures the appearance 
of a different value for the isotope effect would provide support for the hypothesis of two 
competitive mechanisms. 

To investigate this possibility, the mass dependence of the self-diffusion rate for 
zirconium has been measured at four temperatures in the /3-phase. The mass depend- 
ence is related to the isotope effect E at a temperature T by the expression (refs. 11 
and 12) 


E = 




+ (n - l) m 0 \ 1/72 
+ ( n - l) m 0 y 


- 1 


( 1 ) 


where n is the number of ions involved in each diffusion jump, m 
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respective masses of the tracer ions Zr and Zr , 


89 


and mgg are the 


and mg is the mean mass of the 


solvent ion. In the self-diffusion isotope effect measurement the tracers are of the 
same chemical species as the solvent ions, and the correlation factor for diffusion f c 
is related to E by the expression (refs. 13 and 14). 


E = AKf (0 < AK< 1) 
c 


( 2 ) 


where the constant AK represents the relation between the ’'effective" masses of the 
tracer ions required for the derivation of f and the real masses used in the definition 
of E. In general, AK is expected to be close to unity, indicating that the diffusing 
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atom possesses most of the translational kinetic energy during its motion across the 
saddle point of the diffusion barrier. This assumption has been verified experimentally 
for self -diffusion in the metals palladium (face centered cubic) (ref. 15) and iron (face 
centered cubic and body centered cubic) (ref. 10). 

A description of the influence of parent and daughter decay on tracer diffusion pro- 
files by Donald R. Behrendt is given in the appendix. 


EXPERIMENTAL TECHNIQUE 

Cylinders of three pass zone refined zirconium, purchased from the Materials 

Research Corporation, were annealed at 1200 C for 4 hours, under a vacuum of 10" 

torr, in an aluminum oxide tube resistance furnace in the presence of a zirconium foil 

getter. This produced a stable grain size of 5 to 7 millimeters. The ends of each 

cylinder were ground and lightly etched with a solution of 50 percent water, 47 percent 

nitric acid, and 3 percent hydrofluoric acid. An ammonium oxalate solution containing 

89 95 

approximately equal quantities of Zr and Zr was placed on one end of each sample 

89 95 95 

and evaporated to dryness. The Zr and Zr -Nb mixture were obtained from the 

95 

Nuclear Science and Engineering Corporation and an extraction of Zr was performed 

immediately prior to the diffusion anneal using a Dow ex ion exchange column (ref. 16). 

95 95 

The Zr contained less than 1 percent Nb after this procedure. The extraction effi- 
ciency was determined by following the subsequent appearance of the y peak for the 
95 

daughter Nb using a solid state detector with a resolution of 2 keV. 

A pair of cylinders, with their diffusion faces separated by a piece of zirconium 
foil, were wrapped in the same foil and placed inside tantalum cans sealed by electron 
beam welding under a vacuum of 10 ® torr. The diffusion heat treatments at lower tem- 
perature were performed in the furnace already described, and the runs at 1784° and 
1840° C were performed in a tungsten basket furnace. Both furnaces maintained a 
vacuum below 10 -D torr during the diffusion aneals. Temperature was measured to 
within ±0. 1 percent using either platinum - platinum- rhenium or tungsten - tungsten- 
rhenium thermocouples in good thermal contact with the tantalum can. 

After diffusion had occurred, the samples were removed from the tantalum cans and 
each was reduced in diameter (to avoid errors introduced by vapor or surface diffusion) 
and sectioned by lathe or microtome (ref. 17). The individual sections were weighed, 
dissolved in hydrochloric acid, and placed in thin-walled plastic vials. A 3-inch- 

diameter sodium iodide "well” scintillation detector with a resolution of 8 percent was 

89 

used to determine the Zr activity of each section. The y peak at 0. 91 MeV was iso- 
lated by a single channel analyzer, and sufficient counts were recorded to give an ac- 
curacy of 0. 1 percent. 



Figure 1 shows that the y resolution of the detector was sufficient to exclude any 
95 95 89 

significant level of Zr -Nb activity in the Zr channel. 

Samples were counted in random order, and the section having the highest activity 
was used at intervals during the count to determine any changes in the measuring equip- 
ment. The dead time of the system was determined before and after each counting 

session, and correction was made for lost counts, background, and decay during the 

89 

counting. The sections were then stored for 10 half-lives of the Zr (33 days) and 

95 95 

counted for their combined Zr -Nb activity. Correction for the decay during counting 
95 95 

of the Zr -Nb was made by using a decay coefficient for the mixture at the time of 

89 95 95 

counting. Counting to standard solutions containing an identical ratio of Zr /Zr -Nb 

89 95 95 

but having differing total activities showed that measured ratios of Zr /Zr -Nb 
were independent of count rate. 

At the highest temperatures the quantity of daughter produced in the time interval 

95 

between purification of the Zr and the completion of the heat treatment was sufficiently 
small to make correction for the daughter concentration formed during the anneal un- 


Energy, Zr^-’-Nb 9 ' 5 



Figure 1. - Energy spectra for Zr 89 and Zr^-Nb 9 ^ 
mixture using 3- by 3-inch (7.62- by 7.62-cm) 
Nal well scintillation detector. 


o Experimental data determined from 

counts of zr 89 and zr 95 -Nb 95 mixture 
a Data obtained by correction derived 

in appendix 

Theoretical curve for body centered 

cubic metal with simple vacancy exchange 


Abscissa 

scale 



Figure 2. - Isotope effect for zr 89 and zr^-Nb" in 
p-zirconium. 
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necessary. At the lowest temperatures, however, a correction was applied to the 
95 95 

final Zr -Nb count. This was based on a calculation by Behrendt and is given 

95 95 

in the appendix. It is the relative concentration of Zr plus Zr which has decayed 
95 89 

into Nb compared to that of Zr which is required for the isotope effect. Thus the 
95 

presence of Nb daughter does not of itself introduce an error into the results. An 

95 95 

error arises only from the fraction of Nb present which diffused as Nb during the 
diffusion anneal. 

oq qk 95 

Each section was counted twice for both Zr oa and Zr -Nb , and this resulted in 
four plots of 


Inf— against 1 1 - 89 
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V 


d 95 >/ 
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where Ngg is the corrected count for Zr , Ngg represents the count on both Zr and 
Nb 98 , and D is the diffusion coefficient for the appropriate tracer isotope. The pene- 
tration depth and time for the diffusion anneal are given by x and t, respectively. The 
isotope effect data were also plotted as 


N 


In 
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89 

* 

95 
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against (l — )ln (N R q) 
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95 / 
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There was no significant difference in the slopes determined by the two methods, and 
the former expression was used in figure 2. 


CHEMICAL AND RADIOCHEMICAL IMPURITY ANALYSIS 

A mass spectrographic analysis of the original material is given in table I and is 

compared with that determined after a grain growth and high-temperature diffusion 

anneal had been performed. This analysis of the original material differs from the 

manufacturer’s batch analysis, which suggested an overall purity of 99. 99 percent. No 

89 95 95 

radioactive impurities were determined in either the Zr or Zr -Nb y spectra re- 
corded by the solid state detector after the preparation of the diffusion samples. The de- 
tection efficiency was at least 1 percent close to the y peaks and considerable better at 
other y energies. Gamma spectra taken from the first few sections after the diffusion 
heat treatment showed no evidence of slow diffusing impurities, but the lower activity of 
the deeper sections did not permit an exhaustive search for fast diffusing impurities. 


5 



TABLE I. - SPECTROGRAPHIC ANALYSIS OF NIOBIUM AND ZIRCONIUM 


[Results of analysis are expressed as ppm by weight; other elements when present did not 
exceed 5 ppm. ] 


Element 

Isotope 

Analysis of original 
material (Zr) 

effect material 

Analysis after grain growth 
and diffusion (Zr and Nb) 

(a) 

High temperature anneal 
experiment material - analysis 
of original material (Zr) 

C 

35 

120 

97 

N 

130 

36 

66 

O 

140 

1100 

722 

F 

15 

27 

--- 

S 

40 

65 

--- 

Cl 

30 

25 

--- 

Cr 

66 

39 

10 

Fe 

460 

350 

180 

Ni 

4 

8 

--- 

A1 

51 

44 

10 

Hf 

37 

52 

100 

Ta 

2 

4 

— 


a These specimens were given grain growth heat treatment at 1200° C for 4 hr 

under dynamic vacuum of 10 ® torr followed by diffusion heat treatment for 1 hr at 
1784 C in tantalum can sealed under vacuum of 10 torr. 


EXPERIMENTAL RESULTS 

89 

The measured diffusion coefficients for Zr and those of Federer and Lundy 
(ref. 6) for Zr 95 -Nb 95 are shown in figure 3. The diffusion profiles and isotope effect 
plots are given in figures 4 and 2, respectively. The penetration profiles obey the 
following theoretical relation for diffusion from a thin depletable source into a semi- 
infinite solid within the experimental error of the data points: 


C(x, t) 



The initial quantity of tracer is represented by Sq, and C(x, t) is the concentration of 
tracer at point x at time t. The diffusion coefficients calculated from the slopes of 
the penetration profiles are given in table II. 

The isotope effect plots are straight lines for the measurements at the higher tem- 
peratures, but at 942° C there is an increase in the ratio of Zr 89 /Zr 95 -Nb 95 close to 
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a Zr 95 -Nb 95 Federer and Lundy (ref. 6) 1 Diffusion from 

a Nb" Federer and Lundy (ref. 6) r dried deposit of 

o Zr 89 Present report J oxalic acid 



1800 1500 1300 1100 900 

Temperature, °C 

Figure 3. - Plot of log D against temperature for 
Zr 89 , Zr 95 , and Nb 95 diffusing in p-zirconium. 



Figure 4. - Penetration profiles for zr 89 diffusing in 
p-zirconium. 
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TABLE n. - DIFFUSION COEFFICIENTS 


FOR 0- ZIRCONIUM 


Tracer 

Temperature, Time, Diffusion coefficient, 

°C sec cm^/sec 

^ 89 
Zr 

942±1 1.092X10 4 9.88X10" 10 

9.23 

N 

CO 

CD 



1186±1 1.656X10 4 5.25X10" 9 

5.45 

Zr 89 

1784±2 3.630X10 3 2.08X10' 7 

2.21 

Zr 89 

181 5±2 3.510X10 3 2.45X10' 7 

2. 48 

Nb 95 

4 a -10 

1025±20 1.346x10 7. 12x10 

b 2. 03 


a Annealed at 1750° C for 8 hr and diffused after 
30 min at 750° C in a-phase. 

^Annealed at 1750° C for 8 hr and diffused without 
entering a-phase. 


TABLE HI. - SELF-DIFFUSION ISOTOPE FOR /3- ZIRCONIUM 


Temperature, 

°C 

/Dog \ /Dgo \ 

Values of| 1 for each sample Mean value of 11 Temperature 

\ D 95 ) \ D 95 / ratio ’ 

(a) 

942 

-0.002, -0.006, 0.000, -0.003 -0. 001±0. 003 -0. 03±0. 07 

-0.001, 0.001, -0.002, 0.002 

1186 

0.000, -0.002, 0.001, 0.000 0. 000±0.001 0.00±0. 03 

0.001, 0.001, 0.000, -0.001 

1784 

0.001, 0.010, 0.002. 0.000 0. 002±0. 005 0. 06±0. 14 

0.000, 0.005, 0.001, 0.000 

1815 

0.000, -0.001, -0.001, -0.001 0. 000±0. 002 0. 00±0. 06 

0.002, -0.001, -0.001, -0.001 


a Determined by assuming n = 1. 



the specimen surface. It is suspected that the oxalate complex on the surface contri- 
buted impurities which remained in the diffusion zone throughout the diffusion anneal. 

The breakdown of this impurity complex liberated oxygen, nitrogen, and hydrogen in 
quite high local concentrations. However, the very high diffusivities of nitrogen and 
hydrogen (ref. 18) would lead to a rapid reduction in their concentration in a time short 
compared to the tracer diffusion anneal at all temperatures. For oxygen (ref. 18) this 
is true at high temperatures; however, at 942° C the diffusivity of oxygen is a factor 10 
slower than that for zirconium self-diffusion. Thus a high local concentration >1. 0 per- 
cent could explain the anomally in the isotope effect in a region corresponding to the 

diffusion zone for oxygen. The relative diffusion rate D QQ /D Q _ was determined by 

O Cr* o 

least squares analysis for all data points at 1120 , 1784 , and 1840 C for those at 

_5 Q 

greater depths than 4x10 centimeter squared at 942 C. 

From table HI we observe that the isotope effect is less than 0. 2 suggesting that the 

self-diffusion tracers diffuse at a rate which is effectively independent of their mass. 

This behavior is not in agreement with the recent work by the author (ref. 10) on the 

self- diffusion isotope effect in body centered cubic iron, where the mass dependence of 
55 59 

tracers Fe and Fe is consistent with diffusion via a single vacancy exchange mech- 
anism giving an isotope effect of 0.7±0. 06. Before discussing the isotope effect, the 
sources of possible measurement errors, which might give rise to such a result, will be 
reviewed briefly. 

The simple decay technique used to determine the activity distribution of the indi- 
vidual tracers does not raise any problems in the interpretation of the experimental re- 
sults, and our discussion will be limited to an examination of the possible presence and 
influence of impurities during the diffusion anneal. 

Examination of the tracer spectra, both before and after the diffusion heat treat- 
ment, using a solid state detector of 2 keV resolution, showed no characteristic y peaks 

95 

which could be attributed to any impurity other than Nb . Furthermore, serious con- 
tamination by radioactive impurities would be shown by nonlinearity in both the penetra- 
tion profile and isotope effect plots, and it can be seen from figures 4 and 2, with the 

exceptions of the measurements at 942° C, that this did not occur. 

95 95 

The y spectrum of Nb is identical to that of Zr within the 8 percent resolution 

of the sodium iodide counter and the diffusion of a significant quantity of this material 

95 

during the heat treatment would have increased the measured isotope effect, as Nb 

95 95 

diffuses more slowly than Zr . In order to minimize the quantity of Nb present at 

the beginning of the diffusion anneal, a maximum of 1 hour was permitted to elapse be- 

9 5 

tween the purification of the Zr and the start of the diffusion anneal. The quantity of 

95 95 95 

Nb produced before and during the diffusion run can be computed from the Zr -Nb 

transformation curve in the appendix (fig. 5), and it can be shown that correction of the 
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Zr 9 -’ 

Nb 95 



Time, hr 


0 20 40 60 

Time, days 

Figure 5. - Activity of Nb 9 ^and Zr 95 as function of time 
with concentration of Nb 95 zero at time zero. 

95 95 

Zr penetration profile for Nb is insignificant for the high-temperature runs and 
reaches a magnitude of only a few percent at the lowest temperatures. 

An analysis of nonradioactive impurities present in the zirconium before and after 
a high-temperature diffusion anneal shows no significant change in the levels of distri- 
buted substitutional or intersititial impurities with the singular exception of oxygen, 
which is increased by a factor of 8. Work by Peart (ref. 19) on solute diffusion in 
/3-titanium has shown that there is no significant difference in the diffusion properties of 
iodide titanium as compared to commercial 120 titanium. In the former material in- 
dividual substitutional impurities are at levels of a few ppm, while oxygen is at a level 
of around 200 ppm. For the latter material substitutional impurities such as iron, reach 
a level of 2x10 ppm; oxygen is normally 500 ppm and has been doped to levels of 
10^ ppm without any significant change in the measured diffusion coefficient (unpublished 
data obtained at Lewis). Unfortunately, the analysis does not evaluate the impurities 
present in the diffused region during the heat treatment, but it is clear from the above 
work on j3-titanium that quite high levels of contamination would not seriously affect the 
measured diffusion coefficients. 
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The isotope effect measurement at 942° C in this report suggests that a high con- 
centration of interstitial oxygen can alter the isotope ratio where no change in the pene- 
tration profile is detectable but the ratio appears to be altered in a highly localized area 
where the impurity is concentrated. Clearly, when the concentration of oxygen reaches 
a level where the configuration of the metal lattice is significantly changed, the isotope 
effect will not be representative of the pure zirconium lattice, but it seems justifiable to 
assume that, at or below 1000 ppm oxygen and where the pure zirconium lattice is mod- 
erately intact, the relative diffusion rate of two isotopes differing only in their mass 
from each other and the lattice atoms will not be strongly dependent on the presence of 
oxygen interstitials. 

The method of depositing the isotope on the specimen surface was dictated by the 
95 95 

necessity for a rapid Zr , Nb separation and the minimal loss of time before starting 

9 5 

the diffusion anneal. In a diffusion study of Nb in /3-zirconium Federer and Lundy 
(ref. 6) used a similar method for depositing the tracer. The author has made measure- 
ments on the same system, which are described later in this report, using a vacuum 
evaporation technique for depositing the tracer and has found no significant difference in 
the diffusivities measured by the two methods. The presence of barriers to diffusion at 
the surface might be expected in solute diffusion experiments, but such problems do not 
appear to arise in the presence of the oxalate complex and are even less significant in 
self-diffusion experiments where isotopic exchange would guarantee unimpeded entry of 
the tracer into the host lattice. 

From these considerations it is concluded that the results obtained for the isotope 
effect are not artifacts of the experimental method but represent the real mass depend- 
ence of self- diffusion tracers in this system. 


DISCUSSION 

A small value for the isotope effect has been observed for iron diffusing in 
/3-titanium (ref. 20). This result is more difficult to interpret than the self- diffusion 
measurement, but it seems reasonably certain that the major mechanism of diffusion in 
these metals is not simple vacancy exchange. The two mechanism hypothesis has not 
been justified by the results, though it could be argued that the metal melts before nor- 
mal vacancy exchange can successfully compete with the low temperature migration 
process. 

Some years ago an experiment was performed by Kidson and McGurn (ref. 21) in 
which zirconium diffusion samples were annealed for up to 1 hour at 1200° C followed 
by deposition of a tracer without lowering the sample temperature through the at/fi 
phase transition. The diffusion rates obtained by this method were found to be no dif- 
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ferent from those determined by a conventional experiment, and it was assumed that any 
defects introduced at the phase change must anneal out in a time which is small com- 
pared to the length of a typical diffusion anneal. 

In view of the unusual result for the isotope effect which might indicate diffusion via 
short circuit paths such as dislocations, the Kidson and McGurn experiment was re- 
peated using a much higher annealing temperature. 

95 

Two samples of Zr were prepared for deposition of the tracer isotope Nb but were 
first raised to a temperature of 1750° C for 8 hours in an induction furnace under a 
vacuum of 10”® torr. The temperature was then reduced to 1025° ±20° C and held con- 
stant within ±5° C. One specimen was then lowered from the induction coil and cooled to 
750° C (120° C below the transformation temperature). After 30 minutes at this tem- 
perature the specimen was raised to its original position in the heater coil, and the end 

95 

surfaces of both specimens were covered by thin evaporated deposits of Nb . The 
specimens were heated for a further 4 hours and then sectioned by microtome. The dif- 
fusion coefficients measured from the penetration profiles differed by a factor of 3. The 
penetration profiles for these two samples are shown in figure 6. The higher coefficient 

was determined from the specimen which had been passed through the of//3 phase tran- 

89 

sition immediately before the diffusion heat treatment. Zr was temporarily unavail- 
able at the time the experiments were performed, and the solute tracer most closely 
resembling the solvent lattice ion was used. 

The host material was of a lower level of purity than that used for the isotope effect 
experiment, but as enclosure in tantalum cans was not possible, there seemed to be no 



Figure 6. - Penetration profiles for Nb^ diffusing in 
p-zirconium at 1025° C. 
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significant advantage in the use of triple zone refined material. Both diffusion samples 
had almost identical thermal histories, and it would be unreasonable to attempt to ex- 
plain the difference in diffusion properties in terms of differing levels of impurities in 
the two samples. 

It is concluded from this experiment that the a/13 transition in zirconium produces 
some form of short circuit path or complex which is capable of enhancing the solute dif- 
fusion rate and can only be removed by relatively long anneals at temperatures close to 
the melting point. It should be noted that the time required at high temperature is long 
compared to the length of a typical diffusion anneal at these temperatures. 

It is clear that a full understanding of the isotope effect results must await a de- 
tailed examination of this new short circuit phenomenon in zirconium. However, it is 
appropriate to mention some supporting evidence for the presence of a high dislocation 
density in the /3-phase of zirconium extending to much higher temperatures than has been 
proposed on the basis of the observation of time dependent diffusion coefficients for im- 
purity diffusion in /3-titanium close to the a/[ 3 transformation reported by Gregory and 
Askill (ref. 22). 

The suggestion that short circuit paths, dislocations, for example, might be re- 
sponsible for the unusual diffusion behavior of the anomalous body centered cubic metals 
is not new, but there has been considerable difficulty in accepting the large estimates 
of the dislocation density required to explain the observed diffusion coefficients. The ab- 
sence of detectable grain boundary diffusion by autoradiographs even at the lowest tem- 
peratures in the /3-phase of titanium (ref. 19) suggests that the rate of migration along 
grain boundaries and, by implication, also along dislocations does not differ greatly from 
that through the crystal volume. Thus, in order to explain observed diffusion rates at low 

temperature entirely in terms of dislocation diffusion, it is necessary to postulate dis- 

9 11 

location densities in the range 10 to 10 per square centimeter (ref. 23). 

To the author's knowledge there are no direct determinations of dislocation density 
in the |3-phase of either titanium or zirconium, and the only alternative is to examine 
some of the properties of these metals in order to see if any circumstantial evidence 
exists for a high density of dislocations. 

The metals zirconium and titanium both undergo self-deformation at high tempera- 
tures. This appears to involve a gross movement of single crystal regions with respect 
to one another. This phenomenon has only been observed at temperatures above 1600° C 
and has often led to rejection of diffusion samples as it has been impossible to section 
parallel to the original interface because of the relatively small surface area of the dis- 
oriented single crystals. This massive movement of large regions with respect to 
one another could readily absorb large numbers of dislocations and lead to a significant 
decrease in the dislocation density. This behavior could be explained by postulating 
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large dislocation densities immediately after the «//3 transformation which are stable 
until a relatively high temperature is reached. 

Only one electron microscope examination of zirconium is known to the author. 

This work on a -zirconium by Howe, Whitton, and McGurn (ref. 24) suggests that den- 

Q 

sities of the order of 2x10 per square centimeter are typical of this material, but it is 
necessary to show that this dislocation density remains unchanged by passage through 
the phase transformation. 

Glen and Pugh (ref. 25) showed that for a transition from a to ^ to a a single 
crystal of a -zirconium might be expected to transform into one of 57 different a orien- 
tations. Laue photographs show that such a series of transformations invariably re- 
produces the original a orientation. They suggested that dislocations and stacking 
faults were probably responsible for the "memory" associated with the diffusionless 
shear transformation observed in this metal. 

The /3- dislocation densities appear to be stable in the transformation, and high den- 
sities of dislocations may be transferred from one phase to the other. There is clearly 
a necessity for further experimental evidence before short circuit paths can be con- 
sidered an explanation of the diffusion anomaly in /3-zirconium. It is unlikely that it will 
be possible to obtain electron microscope transmission pictures of dislocations in the 
/3-phase, but work is at present being undertaken on a detailed investigation of the 
structure of the a- and /3-phases of zirconium using X-ray techniques (unpublished data 
from G. V. Kidson). An experimental examination of the influence of high temperature 
anneals on self-diffusion is also underway at the NASA Lewis Research Center. It is 
hoped that an investigation of the kinetics of high-temperature annealing will make it 
possible to determine the real self- diffusion temperature dependence in the absence of 
large concentrations of short circuit defects. 

The isotope effect for diffusion in a dislocation matrix of the kind suggested is dif- 
ficult to estimate. If the problem is considered as diffusion along an idealized edge 
dislocation, which might be the case for the strong segregation of an impurity tracer, 
then we should expect an isotope effect corresponding to the inverse square root of the 
mass. This, however, seems to be an unrealistic assumption as it does not consider 
the disorder existing in a real dislocation, and it might be expected that the opportunity 
for relaxation close to the dislocation core would allow considerable kinetic energy 
transfer between the diffusing ions and those ions forming the diffusion saddle. This 
could give a value of AK < 1. Alternatively, considerable cooperative motion of ions 
may occur during each diffusion exchange, leading to a dilution of the mass dependence 
of the tracer ions. At least seven additional ions would be required to account for the 
observed isotope effect. Either factor or a combination of both could provide an ex- 
planation for the small isotope effect. 
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CONCLUDING REMARKS 


89 95 

The isotope effect has been determined for Zr and Zr diffusing in /3-zirconium. 
At the temperatures 942°, 1186°, 1784°, and 1815° C the isotope effect was determined 
to be less than 0. 2. An annealing experiment to determine the importance of the dif- 
fusionless shear transformation in contributing diffusion short circuit paths to the 
/3-phase of zirconium has indicated that a threefold reduction in the measured diffusion 
coefficient at 1025° C can be achieved by a heat treatment of several hours close to the 
melting temperature. 

The unusual Arrhenius relation for zirconium, titanium, and hafnium and the small 
isotope effect are tentatively attributed to diffusion on a short circuit network. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 8, 1970, 

129-03. 
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APPENDIX - INFLUENCE OF PARENT AND DAUGHTER DECAY 
ON TRACER DIFFUSION PROFILES 
by Donald R. Behrendt 

The kinetics of parent- daughter decay are given by the expressions 

-X t 

N p = N 0 e P 



where Nq is the initial disintegration rate for the pure parent isotope at time t = 0, 

Np is the disintegration rate for the parent at time t, and is that for the daughter. 

The symbol X is used to denote the decay constant associated with either the parent 

(X ) or the daughter (X,). Figure 6 shows the time dependence of the count rates for Zr^® 
P qr a n n 

and Nb . It is clear that at 1784 and 1815 C, where the diffusion coefficients for 
95 95 

Nb and Zr are approximately equal, the heat treatments were completed in such a 

95 95 95 

short time that no significant correction to the Zr -Nb count is required for Nb 
formed prior to or during the diffusion anneal. 

In order to calculate the influence of daughter diffusion during the longer heat treat- 
ments required at 942° and 1186° C, an equation for the penetration profile of the com- 
bined parent-daughter combination is derived as follows. Consider a one -dimensional 
diffusion geometry where n p (x, t) and n^(x, t) are, respectively, the volume concentra- 
tions of the parent (the tracer) and the daughter decay product as a function of the dis- 
tance x, and the time t, where x is measured from the surface of a semi- infinite 
sample. The dependent variables will then be independent of y and z. 

Let the parent have a decay constant X p and a diffusion coefficient D p . Similarly, 
the daughter has a decay constant X^ and a diffusion coefficient D^. The concentration 
of the parent is then a solution to the equation 


0n (x, t) 9 n (x, t) 

_P = D 2 X n (x, t) 

0t P o P P 

^ 9x Z 


(Al) 


16 


. 



!! 

L 

t 


This is the usual one-dimensional diffusion equation except for the last term, which 
takes into account the rate of decay of the parent. If the parent has an initial uniform 
concentration n n in a region of thickness a at the surface of the sample, the bound- 
ary conditions are 


n p (x, 0) = n pj 0 0 < x < a 


x > a 


3n (°, t) 

— 2- = o lim n (x, t) = 0 

3x p 


x — °° 


(A 2) 


Similarly, the differential equation for the concentration of the daughter is 


0n d (x, t) 3‘ i n d (x, t) 


= D 


at 


+ x n n J x > t) - x d n d (x ’ 


3x 


P P 


(A3) 


and the boundary conditions are 


n d< x >°> = n d,0 0£ 


x < a ^ 


= 0 


x > a 


3n d (0, t) 


V 


3x 


= 0 lim n d (x, t) = 0 


x — 




(A 4) 


where n d Q is the initial uniform concentration of the daughter in the thickness a at 
the surface 

If the thickness, a, of the initial radioactive plating is small compared to the pene- 
tration depth, the solutions to these equations are 
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(A6) 


where N and N, are, respectively, the initial surface concentrations of the parent 

r 

and daughter and are related to the initial volume concentrations by 


N p = n P o a ’ lim 


N d ~ n dO a > lim 



(A7) 


The total counting rate r^ of a thin section of thickness AZ and area A including 
both parent and daughter decays is 


r pd = ^ AA '[v P V x ’ t) 


+ f 


d X d n d( x ’ t )] 


(A8) 


where J/ represents the efficiency of the counter for counting the gammas and f and 
f^ are the respective fractions of the decays of the parent and daughter which result in 
the measured gamma rays. 

-7 

As an example, the concentration n^(x, t) was calculated for equal to 1x10 
reciprocal seconds, and t equal to 2x10^ seconds for several values of Dp and D^. 

Figure 7 is a plot of the log of fp^/i/FA A Z as a function of x^/t and shows these 
results. The short-dashed curves in these figures represent the usual diffusion profile 
obtained from the counting rate of the parent only. Each of the other curves then repre- 
sents the diffusion profile obtained by summing the counting rates of the parent and 
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Figure 7. - Reduced counting rate as function of diffusion profile parameter. 


daughter. The curve for D p = 1x10 ® square centimeter per second and D d = 8x10*° 
square centimeter per second represents the measured penetration profile for Zr^-Nb^ 
diffused into zirconium at 942° C for 4 hours. 

Figure 7(a), where D p > D d , shows clearly that the slopes derived from the solid 
curves and hence the apparent diffusion coefficients are only a few percent smaller than 
the true value given by the reciprocal slope of the dashed line representing the parent 
alone. However, the influence of the daughter is much more marked in figure 7(b) 
where D d > D . In the worst case presented, diffusion coefficients calculated from a 
straight line drawn through curves with D d /D p of 10 or 100 could overestimate D 
by 25 percent. 


••i 


19 



REFERENCES 


1. Lundy, Ted S. : Diffusion in Body-Centered Cubic Metals Zirconium, Vanadium, 

Niobium, and Tantalum. Rep. ORNL-3617, Oak Ridge National Lab. , June 1964. 

fiO 

2. Askill, J. : The Effect of an Argon Environment on the Diffusion of Co in Molyb- 

denum. Phys. Stat. Sol., vol. 9, 1965, pp. K113-K116. 

3. Kidson, G. V.; and Ross, R. : Radioisotopes in Scientific Research. Vol. 1. 

Pergamon Press, 1958, p. 185. 

4. LeClaire, A. D. : The Application of Diffusion Theory to the Body Cubic Centered 

Structures. Diffusion in Body-Centered Cubic Metals, Proceedings of Gatlinburg 
Conference. American Society for Metals, 1964, pp. 3-26. 

44 48 

5. Murdock, J. F. ; Lundy, T. S. ; and Stansbury, E. E. : Diffusion of Ti and V in 

Titanium, Acta Met. , vol. 12, no. 9, Sept. 1964, pp. 1033-1039. 

95 95 

6. Federer, J. I.; and Lundy, T. S. : Diffusion of Zr and Cb in Boc Zirconium. 

Trans. AIME, vol. 227, no. 3, June 1963, pp. 592-597. 

181 

7. Winslow, F. R. ; and Lundy, T. S. : Diffusion of Hf in Bcc Hafnium. Trans. 

AIME, vol. 233, no. 9, Sept. 1965, pp. 1790-1791. 

8. McLean, Donald: Grain Boundaries in Metals. Oxford University Press, 1957. 

9. Graham, D. : A Summary of the Impurity Diffusion to Beta Phase of Titanium. Dif- 

fusion in Body Centered Cubic Metals, Proceedings of Gatlinburg Conference. 
American Society for Metals, 1964, pp. 27-35. 

10. Graham, D. : Mass Dependence of Self-Diffusion in Iron. J. Appl. Phys. , vol. 40, 

no. 6, May 1969, pp. 2386-2390. 

11. Schoen, A. H. : Correlation and the Isotope Effect for Diffusion in Crystalline 

Solids. Phys. Rev. Letters, vol. 1, no. 4, Aug. 15, 1958, pp. 138-140. 

12. Vineyard, George H. : Frequency Factors and Isotope Effects in Solid State Rate 

Processes. J. Phys. Chem. Solids, vol. 3, no. 1/2, 1957, pp. 121-127. 

13. Barr, L. W. ; and LeClaire, A. D. : Correlation Effects in Diffusion. Proc. Brit. 

Ceram. Soc. , vol. 1, 1964, pp. 109-128. 

14. Mullen, James G. : Isotope Effect in Intermetallic Diffusion. Phys. Rev., vol. 121, 

no. 6, Mar. 15, 1961, pp. 1649-1658. 

15. Peterson, N. L. : Isotope Effect in Self-Diffusion in Palladium. Phys. Rev. , 

vol. A136, no. 2, Oct. 19, 1964, pp. 568-574. 


20 



16. Huffman, E. H. ; Iddings, G. M. ; and Lilly, R. C. : Anion Exchange of Zirconium, 

Hafnium, Niobium and Tantalum in Hydrochloric Acid Solutions. J. Am. Chem. 
Soc., vol. 73, no. 9, Sept. 1951, pp. 4474-4475. 

17. Graham, D. : Tracer Diffusion Studies on Metals Using a Microtome. Rev. Sci. 

Instr. , vol. 40, no. 7, July 1969, pp. 897-901. 

18. Miller, George L. : Zirconium. Academic Press, 1957. 

19. Peart, Reginald F. : Diffusion in Body Cubic Centered Titanium and Titanium 

Alloys. PhD Thesis, Reading University, 1959. 

20. Gibbs, G. B. ; Graham, D. ; and Tomlin, D. H. : Diffusion in Titanium and 

Titanium -Niobium Alloys. Phil. Mag., vol. 8, no. 92, Aug. 1963, pp. 1269-1282. 

21. Kidson, G. ; and McGurn, J. ; Self-Diffusion in Body-Centered Cubic Zirconium. 

Can. J. Phys., vol. 39, no. 8, Aug. 1961, pp. 1146-1157. 

22. Gregory, A. G.; and Askill J. : Time Dependence and Other Effects on the Tracer 

113 

Diffusion of Sn in Titanium. Phil. Mag., vol. 12, no. 119, Nov. 1965, 
pp. 901-906. 

23. Graham, Donald: Atomic Diffusion in the Titanium-Niobium System. PhD Thesis, 

Reading University, 1963. 

24. Howe, L. M. ; Whitton, J. L. ; and McGurn, J. F. : Observation of Dislocation 

Movement and Interaction in Zirconium by Transmission Electron Microscopy. 
Acta Met. , vol. 10, no. 9, Sept. 1962, pp. 773-787. 

25. Glen, J. W. ; and Pugh, S. F. : The Effect of Phase Transformations on the Orien- 

tation of Zirconium Crystals. Acta Met. , vol. 2, no. 3, May 1954, pp. 520-529. 


NASA-Langley, 1970 26 E-5094 


21 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


official business 


FIRST CLASS MAIL 


POSTAGE and fees paid 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


04U 001 A2 51 3DS 70195 00903 

AIR FORCE WEAPONS LABORATORY /WLOL/ 
KIRTLAND AFB, NEW MEXICO 87117 


ATT E. LOU BOWMAN, CHIEF, TECH. LIBRARY 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


■ ■','The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of, ..phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.’’ 

— National Aeronautics and Space Act of 1958 

1 4 

NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 


published in a foreign language considered 
to merit NASA distribution in English. 


SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 


TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



